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This paper theoretically discusses the spin current in spin-triplet superconductor / insulator 
/ spin-triplet superconductor junctions. At low temperatures, a midgap Andreev resonant state 
anomalously enhances not only the charge current but also the spin current. The coupling between 
the Cooper pairs and the electromagnetic fields leads to the Frounhofer pattern in the direct current 
spin flow in magnetic fields and the alternative spin current under applied bias-voltages. 
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Although the supercurrent usually refers to the dis- 
sipationless charge flow in superconductors or the mass 
flow in superfluid He, the spin flow carried by spin-triplet 
Cooper pairs is also undoubtedly a supercurrenti. These 
supercurrents have a common feature; the spatial gradi- 
ent of the order parameter drives the supercurrent. It is 
well known that the charge and mass supercurrents are 
possible under the spatial gradient in the macroscopic 
phase ip. On the other hand, the spatial gradient of d- 
vector causes the spin supercurrent. In the mean-field 
theory of superconductivity, d characterizes the order pa- 
rameter of spin-triplet pairs as A(r ) = id(f ) • ero^e^ 
with <7j for j = 1 — 3 are the Pauli's matrices. In 3 He, 
the angular momentum vector of a Cooper pair points the 
direction normal to the surrounding wall. The spin cur- 
rent is expected near the curved wall because d and the 
momentum vector align to each other due to the dipole- 
dipole interaction. Through the spin dynamics in super- 
fluids, the spin current is detected by nuclear magnetic 
resonance experiments^^. 

In bulk superconductors, however, the spin current is 
usually not expected because the weak spin anisotropy 
fixes d homogeneously in a certain direction of the crys- 
tal lattice. So far the generation of the spin flow has 
been discussed in curved structures of superconductors^ 
and superconducting weak links^. To realize the spin 
current in experiments, we should consider simpler struc- 
tures such as superconductor / insulator / superconduc- 
tor (SIS) junctions of triplet superconductors. Gener- 
ally speaking, spin-triplet pairs are characterized by odd 
parity p- or /-wave symmetry. In SIS junctions of such 
unconventional superconductors, a midgap Andreev res- 
onant state (MARS) governs the low energy transport 7 . 
So far it is pointed out that the MARS drastically en- 
hances the electric charge transport but it suppresses the 
thermal transport 8 . Effects of the MARS on the spin 
transport is still an open question. In contrast to spin- 
triplet pairs of 3 He, Cooper pairs in superconductors have 
the charge degree of freedom which couples with elec- 
tromagnetic fields. This feature may also enables us to 
switch the spin current by using electric fields and/or 
magnetic fields. 

In this paper, I calculate analytically the spin current 
in SIS junctions based on the mean-field theory of super- 
conductivity. When the MARS is forming at the junction 



interface, the low-temperature anomaly appears in the 
spin current as well as it does in the charge currenlsSiifl. 
The spin current shows the Fraunhofer pattern in applied 
magnetic fields. I also show that the applied bias- voltage 
across junctions causes the alternating current of spin. 

Throughout this paper, we take the unit of h = ks = 
c = I , where ks is the Boltzmann constant and c denotes 
the speed of light. The vectors in the real and momentum 
space are indicated by • — while those in the spin space 
are described by bold-italic characters. 
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FIG. I: (Color online) (a) A schematic figure of a SIS junc- 
tion. The large arrows in superconductors represent d which 
is a unit vector in the spin space. In (b), we propose a SIS 
junction using Sr2Ru04, where d is parallel to the c axis of 
the crystal. 



Electronic states in SIS junctions are described by the 



2 



Bogoliubov-de Gennes (BdG) equation, 
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h(r, f ') = S(r- r ') ^ + K(r ) - /i a , (2) 



and hole branches, respectively. Secondly I calculate the 
Andreev reflection coefficients from the boundary condi- 
tions, 



2z o fc F * fl (0,p), 



where Dp = V? — ieA? with Af? being the vector poten- 
tial, ■ • ■ indicates 2x2 matrix describing spin space, &o is 
the unit matrix and /i is the Fermi energy. The insulat- 
ing barrier at x = is described by Vb5(x). In uniform 
superconductors, the pair potential is given by 
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The vector k — (k,p) represents a wave numbers on the 



Fermi surface (i.e. 



k 2 + p 2 



= hp), where k and p are 



wave number in the direction of the current and that in 
the directions transverse to the current, respectively. I 
consider a similar model as that in Ref. 0. We assume 
that the two superconductors are identical to each other 
except for the directions of the unite vectors dj (j = L 
or R) as shown in Fig.^a), where a is an orientation an- 
gle between d^ and d R . In unconventional superconduc- 
tor junctions, the characteristic behaviors of Josephson 
charge current are very sensitive to dependences of the 
pair potential on wavenumbers. When Ak,p and A—k,p 
have opposite signs to each other, the MARS forming at 
junction interfaces causes the anomalous charge trans- 
port at low temperature o 7 ' 12 . In this paper, I consider 
the two typical situations: A_&^ = vAk,p with v = ±1. 
The MARS forms (does not form) for v = — 1 (u = 1). 
Within p-wave symmetries, v — 1 (v — —1) corresponds 
to p y - (p x -) wave symmetry. At first I solve the BdG 
equation in the two superconductors independently and 
obtain the wave functions 
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where u(v) = y/(oj n + (-)n)/(2w n ), = y/u* + |A £ | 2 , 
u = ua , v = uA?/|Ag|, and ui n = (2n+ 1)T is the Mat- 
subara frequency. The incident amplitudes of a quasipar- 
ticle in the electron (hole) branch are denoted by diagonal 
2x2 matrix a (b). In * L A and B (C and D) 

represent the amplitude of outgoing wave in the electron 



with zq = Vb /{kp /to). The Andreev reflection coeffi- 
cients are defied by the off-diagonal elements of the ma- 
trix relation below 
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The calculated results of the Andreev reflection coeffi- 
cients are summarized as 
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K = [(a c — aiCi sin 2 a)<7 + (a ci — Coai)n-&] /hih 2 , 
a o =k| 2 /-(^3+cosa - /+) + \t\ 2 (f + vcosa - g+)hi, 
a a =ivj_ [\r\ 2 g-f- + |t| 2 /ii] , 
c o =h\ {7 — (1 + v) cos a} 

— |r| 2 ft-2 {g+ — u 2 v 2 j — u 2 v 2 (l + v) cos a} , 
ci =-i(l-i/) [hi + |r| 2 /i 2M 2 w 2 ] , 
hi =u 4 e^ + f 4 e-^ 
/i2 =2(cosyj — j/cosa), 
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with /j- = ?i ± i/u , g± 



e - '^^, (p = (p L - <p R , n = d R x d L , and 3 = S + S_. 
The normal transmission and reflection probabilities for a 
channel characterized byp are given by \t\ 2 = k 2 / \z^k 2 F -\- 
k 2 ) and \r\ 2 = z^kp / '{z 2 k 2 F + fc 2 ), respectively. Finally 
I obtain the spin current in SIS junctions based on the 
formula^ 
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The electric Josephson current is also given by <x/2 — » 
— e<To in Eq. Q 13,14 . The spin polarized in the direction 
of n || z flows through the SIS junction as shown in 
Fig. ^a). The electric current J e and the z component 
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of the spin current J* in the x direction result in 
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where Ao is the amplitude of the pair potential at zero 
temperature. The summation in terms of p means the 
summation of the current over propagation channels. It 
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FIG. 2: (Color online) The temperature dependence of the 
spin current for a = 0.1-k and Zo — 10. The results are nor- 
malized by the maximum value of the spin current at zo = 
and T = 0. 

is easy to confirm that J s (J e ) is an odd function of a 
(p). At zo — 0, J e and J* are independent of v be- 
cause the transmission probability is unity (i.e., \t\ 2 = 1). 
In what follows, we consider the low transmission limit 
(i.e., |£| 2 <C 1). We show the maximum amplitude of 
J I as a function of temperatures in Fig. [3 for z$ = 10 
and a = 0.17T. The maximum values are calculated from 
the current-phase (J| — ip) relation and are normalized 
by Jo = Yle l^fel/^ol sin(a/2)| which corresponds to the 
critical current at zq = and T = 0. For v = 1, the tem- 
perature dependence of the spin current is described by 
the Ambegaokar-Baratoff formula 15 and the spin current 
saturates at low temperatures as shown in Fig. EI On the 
other hand for v = — 1, the spin current increases rapidly 
with decreasing temperatures because the MARS forms 
at the interface. The resonant transmission through the 
MARS causes the low-temperature anomaly not only in 
the charge transport^iifii but also in the spin transport. 



This is because the Andreev reflection carries both elec- 
tric charges and spins as shown in Eq. (JJJ. In fact, we 
find the relation J e (a, ip) = Jg(<p, a). Interchange of the 
angle in the gauge space p> and the angle in the spin space 
a connects the two the supercurrents, which implies the 
duality of the charge and spin. 

In contrast to spin-triplet pairs of 3 He, Cooper pairs 
in superconductors couple with the electromagnetic field 
through the charge degree of freedom of an electron. As a 
result, the spin current is expected to be sensitive to elec- 
tromagnetic fields. Here I first consider the spin current 
in SIS junctions under magnetic fields. It is well known 
that the electric Josephson critical current exhibits the 
Frounhofer pattern under magnetic fields. We apply the 
widely accepted wajii& to Eq. © and calculate the criti- 
cal spin current near the critical temperature as 
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where a magnetic field B is parallel to the z direction, 
$ = d n WB with do/ 2 being a penetration depth of 
magnetic fields and $0 is the flux quantum. The crit- 
ical current of spin shows exactly the same Frounhofer 
pattern as that of the charge current. For v = — 1, 
the period of oscillations becomes 2<&o in low temper- 
atures due to the resonant transmission through the 
MARS. Next I discuss the spin current under applied 
bias-voltage. The ac Josephson effect is also well known 
in SIS junctions of spin-singlet superconductors. The 
macroscopic phase obeys the Josephson's equation of mo- 
tion dip/dt — 2eVbias under the bias- voltage T^ias across 
the junction. The spin current under the bias-voltage 
near the critical temperature becomes 



Jg oc cos(2eVbi as £ + Po) sin a. 



(13) 



The bias-voltages generate the alternating spin current 
in SIS junctions. The gauge coupling of the Cooper pair 
results in Eqs. 1|12[> and l|13f) which represent the specific 
properties of the spin current in superconductor junc- 
tions. In addition, it may be possible to control the spin 
current by using these properties. 

To observe the spin current in superconducting mate- 
rials or junctions, we have to overcome two difficulties 
in experiments: the generation and detection of the spin 
current. To generate the spin current, we propose a SIS 
junction of the spin-triplet superconductoi— Sr2Ru04 as 
shown in Fig. ^b). The bending junction is required to 
have the spatial gradient of d. In Sr2Ru04, d homoge- 
neously aligns in the parallel direction to the c axis of 
the crystal lattice. The fabrication of such bending SIS 
junction may be possible because the crystal growth rate 
within the ab plane is much faster than that in the c axis. 
The grain boundary separates the two thin films of single 
crystal. At present, the chiral p x ± ip y -w&ve symmetry 
is the promising candidate of the pairing symmetry in 
Sr2Ru04. The formation of the MARS depends on inci- 
dent angles of a quasiparticle into junction interface. As 
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a result, the MARS is weakly forming at the interface. 
The spin current has an intermediate quality between the 
spin current with v = 1 and that with v = — 1. Although 
Eqs. ©, |Q and l|13fl are derived in the ballistic junc- 
tions, these behaviors persist even in realistic disordered 
junctions. This is because the MARS causes the anoma- 
lous Josephson effect^ in the p-wave symmetry. The 
electromagnetic fields modulate the spin current as dis- 
cussed in Eqs. |(T1|| an d C3- At present, unfortunately, 
there is no good idea of directly observing the spin flow. 
Actually, the spin Hall effect*^ is confirmed by the spin 
accumulation which is a result of the spin flow. The ap- 
plication of the spin-orbit device^ may be possible to 
detect the spin current in superconducting junctions. 

In summary, I have studied the spin supercurrent in 
superconductor / insulator / superconductor junctions, 
where superconductors are in spin-triplet unitary states. 
On the basis of the current formula, I calculate spin and 
charge current analytically from the Andreev reflection 



coefficients of the junctions. The resonant transmission 
through the midgap Andreev resonant state causes the 
low-temperature anomaly of the spin current. The gauge 
coupling of the Cooper pairs with electromagnetic field 
results in the Frounhofer pattern of the spin current un- 
der magnetic fields and the alternating spin flow in the 
presence of bias-voltages. The application of the ob- 
tained results to realistic junctions consisting of Sr2Ru04 
is briefly discussed. 
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